Interference alignment-based sum capacity bounds for 
random dense Gaussian interference networks 



Oliver Johnson* Matthew Aldridge^ Robert Piechocki* 

July 29, 2009 



Abstract 

We consider a dense K user Gaussian interference network formed by paired trans- 
mitters and receivers placed independently at random in a fixed spatial region. 
Under natural conditions on the node position distributions and signal attenua- 
tion, we prove convergence in probability of the average per-user capacity Cs/K to 
log(l+2SNR). The achievability result follows directly from results based on an 
interference alignment scheme presented in recent work of Nazer et al. Our main 
contribution comes through an upper bound, motivated by ideas of 'bottleneck 
capacity' developed in recent work of Jafar. By controlling the physical location 
of transmitter-receiver pairs, we can match a large proportion of these pairs to 
form so-called e-bottleneck links, with consequent control of the sum capacity 



1 Introduction and main result 

1.1 Interference networks and bottleneck states 



Recent work of Jafar [9] made significant progress towards what is referred to as 'the 
holy grail of network information theory', namely the calculation of the capacity of 
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arbitrary Gaussian interference networks. Jafar proves convergence in probability of 
the averaged sum capacity of certain dense Gaussian interference networks. Although 
results contained in the paper [9] made significant progress with this problem, the results 
were described under the constraint that each direct link had the same fading coefficient 
VSNR - a constraint that we relax in this paper. 

In Lemma 1 of [U] , Jafar showed that a two- user Gaussian interference channel with one 
of the cross-link strengths INR = SNR has sum capacity exactly equal to log(l + 2SNR). 
Jafar described such a configuration as an example of a 'bottleneck state', in that altering 
the other cross-link strength does not affect the capacity. Jafar went on to define the 
concept of an e-bottleneck link - that is, a cross-link in a two-user channel with capacity 
within e of log(l + 2SNR). He considered a model of large networks, where each value 
of SNR is fixed, and each INR is sampled independently from a fixed distribution. Jafar 
argues that with probabilty 5, each INR lies in the range such that the corresponding 
two-user channel becomes an e-bottleneck. He uses an argument based on Chebyshev's 
inequality to deduce that 



It is perhaps surprising that the existence of a positive proportion of e-bottleneck links 
implies accurate probabilistic bounds on the sum capacity Cs of the whole network. 
However, we might regard it as analogous to the so-called 'birthday paradox' - although 
each cross-link has a probability 5 of being in an e-bottleneck, there are K(K — 1) cross- 
links that can have this property, so as K tends to infinity, the number of cross-links 
with this property becomes much larger than K. 

In this paper, we show that results such as Equation (PQ) in fact hold more generally, 
in cases where transmitter node positions T 1; . . . , Tk and corresponding receiver node 
positions R±, . . . , Rk are chosen independently at random in a region of space T>. While 
exact expressions for the capacity of networks with arbitrarily placed nodes remain 
elusive, in Theorem 11.51 we prove convergence in probability for the average per-user 
capacity C^/K of such dense network configurations. This may perhaps be analogous 
to the fact that Shannon [14J proved that the average code performed well, while it 
remains a significantly harder problem to establish good performance for a particular 
family of codes. In other words, the randomness we add to the model helps us, rather 
than making things harder. The intuition is that in a dense network of points, a large 
proportion of nodes can be put together pairwise to form e-bottleneck links (see Figure 



While it would be possible to adjust individual transmitters' powers to force each SNR to 
become equal (as assumed by |9]), individual power constraints make this undesirable. 
Further, Jafar assumed that the INR are independent and identicaly distributed, a 




(1) 
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Figure 1: A dense network with K = 6 transmitter-receiver pairs placed on the square 
[0, l] 2 , with e-bottleneck links emphasised. The distances from T 5 to R$, R$ to T\ and 
T\ to R\ are all approximately equal, similarly for T 4 to i? 4 , _R 4 to T 2 and T 2 to R 2 . 
Transmitter-receiver pairs (T 3 , R 3 ) and (T 6 , i? 6 ) are not matched into e-bottleneck links. 



property that would be lost if user powers were scaled in this way. 
1.2 Interference Alignment 

The concept of interference alignment first appeared in [2], and represents a depar- 
ture from the recent paradigm of random (or rather pseudo-random) code construction. 
Pseudo-random codes (e.g. Turbo and LDPC codes) have revolutionised point-to-point 
communications. However, since multi-terminal networks are interference-limited rather 
than noise-limited, unstructured (pseudo-random) codes are not suitable. 

Interference alignment advocates a collaborative solution. Each receiver divides its sig- 
nalling space (space/time/frequency/scale resources) into two parts; one for the signal 
from the intended transmitter, and the second acts as a waste bin. The encoding is 
structured is such a way that the transmitted signal from each of the K transmitters is 
seen in the clear space for the intended receiver, and at the same time, it falls into the 
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waste bin for all other receivers. In such a scenario, the network is no longer interference 
limited. In the first such scheme Cadambe and Jafar [2] showed that spatio/temporal 
beamforming in the SNR — > oo regime allows to achieve for each user pair a "capacity" 
equal to half that of the single user channel. A different idea recently appeared in 
which does not require SNR — > oo. Consider the following channel pair: 



1 -1 
1 1 
-1 1 



H b 



1 -1 
1 1 
-1 1 



(2) 



The channels are constructed in such a way that the simple sum of the received signals 
leads to interference alignment since H a + Hb = 21. In more general channels, such as 
those with Rayleigh fading, one needs to code over sufficiently long time intervals to 
observe and match complementary matrix pairs. 



1.3 Node positioning model 

We believe that our techniques should work in a variety of models for the node positions. 
We outline one very natural scenario here. 

Definition 1.1 Consider a fixed spatial region T> G FP , with two probability distribu- 
tions Ft and Fr supported on T>. Given an integer K , we sample the K transmitter 
node positions Ti, . . . ,Tk independently from the distribution Ft- Similarly, we sample 
the K receiver node Ri, . . . , Rk positions independently from distribution Fr. We refer 
to such a model of node placement as an TID network'. 

Equivalently, we could state that transmitter and receiver positions are distributed ac- 
cording to two independent (non-homogeneous) Poisson processes, conditioned such that 
there are K points of each type in T>. We pair the transmitter and receiver nodes up so 
that Ti wishes to communicate with Ri for each i. We make the following definition: 

Definition 1.2 We say that transmitter and receiver distributions Ft and Fr are 'spa- 
tially separated' if there exist constants C sep < oo and D sep < oo such that for T ~ P T 
and R ~ the Euclidean distance d(T, R) satisfies 

F(d(T, R)<s)< C scp s Dscp for all s. (3) 

We argue in Lemma 12.21 below that a wide range of node distributions Ft and P# have 
this spatial separation property, which allows us to control the tails of the distribution 
of SNR, and hence the maximum value of SNR in Lemma 12.31 
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1.4 Transmission models 



For simplicity, we first describe our results in the context of so-called 'line of sight' com- 
munication models, without multipath interference. That is, we consider a model where 
signal strengths decay deterministically with Euclidean distance d according to some 
monotonically decreasing continuous function f(d). We make the following definition, 
which complements the definition of spatial separation given in Definition 11.21 

Definition 1.3 We say that the signal is 'decaying at rate a' if there exists Cdec < oo 
such that for all d 

f(d) < C dec d~ a . (4) 

Standard physical considerations imply that all signals must be decaying at some rate 
a > D, where D is the dimension of the underlying space. Tse and Viswanath [Tt)J 
Section 2.1] discuss a variety of models under which this condition holds for different 
exponents a. 

We define the full action of the Gaussian interference network: 

Definition 1.4 Fix transmitter node positions \T\, . . . ,7V} G T> and receiver node po- 
sitions {R\, . . . , Rr} G T>, and consider Euclidean distance d and attenuation function 
f. For each i and j, define INR^- = f(d(Ti,Rj)). For emphasis, for each i we write 
SNRj for INRjj . 

We consider the K user Gaussian interference network defined so that transmitter i 
sends a message encoded as a sequence of complex numbers Xj = PQ[1], . . . ,Xi[N]) to 
receiver i, under a power constraint -h ^2 n=1 \Xi[n]\ 2 < 1 for each i. The nth symbol 
received at receiver j is given as 

K 

YM = ^exp(z^M) v ^NR~X 4 [n] + Z>], (5) 
i=i 

where Zj[n] are independent standard complex Gaussians, and 4>ij[n] are independent 
U[0, 2n] random variables independent of all other terms. The INR^ remain fixed, since 
the node positions themselves are fixed, but the phases are fast fading. 

We write Sij for the random variables |log(l + 2INRy), which are functions of the 
distance between % and Rj. In particular, since the nodes are positioned independently 
in Definition 11.11 under this model the random variables Sa = \ log(l + 2SNRj) are 
independent and identically distributed. 

In Section 14.21 we explain how our techniques can be extended to apply to more general 
models, in the presence of random fading terms. 
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1.5 Main result: convergence in probability of Cy,/K 



We now state the main theorem of this paper, which proves convergence in probability 
of the averaged capacity, under the model of node placement described in Definition 11.11 
and the model for signal attenuation described in Definition 11.41 For the sake of clarity, 
we restrict our attention to the case where P# and Pt are uniform, though we discuss 
later to what extent this assumption is necessary. 

We restrict to bounded regions T> with a smooth boundary - the sense of this smoothness 
will be made precise in the proof of Theorem 11.51 Theorem 11.51 will certainly hold for 
squares and balls T> = [0, 1} D and T> = {x : d(x, 0) < 1}, and indeed for any convex and 
bounded polytopes (with finite surface area). Essentially we require that the boundary 
of T> x T> has Hausdorff dimension < 2D — 1, which is very natural. 

Theorem 1.5 Consider a Gaussian interference network formed by K pairs of nodes 
placed in an IID network, with the signal decaying at some rate a > D. If distributions 
Ft and Pr are both uniform on a bounded region T> with smooth boundary, then the 
average per-user capacity C^/K converges in probability to ^Elog(l + 2SNR) ; that is 



lim P 



^-I E log(l + 2SNR) 

K Z 



> e = for alle> 0. 



Proof We break the probability into two terms which we deal with separately. That 
is, writing E = ES U = |Elog(l + 2SNR), 



P 



— — E 
K 



> e )=P(§-£<- e )+p(§-£> e ). (6) 



Bounding the first term of (EJ) corresponds to the achievability part of the proof. Bound- 
ing the second term of (EI) corresponds to the converse part, and represents our major 
contribution. 

The first term of can be bounded relatively simply, using an achievability argument 
based on an interference alignment scheme presented by Nazer, Gastpar, Jafar and 
Vishwanath [11]. Theorem 3 of [Hj implies that the rates R[i] = 1/2 log(l+2SNRj) = Sa 
are achievable. This implies that Cs > Ylf=i This allows us to bound the first term 
in Equation (jSJ) as 

F^-E<-e^j < ¥ ( IM^3 <_ e j, (7) 

Note that Lemma 12.21 below implies that the spatial separation condition holds in the 
setting of Theorem 11.51 Hence the conditions of Lemma [2.31 hold, implying that Su has 
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finite variance by Equation ffT2"]) . This means that Equation ((7|) can be bounded by 
Var (So) /(Ke 2 ), and tends to zero at rate 0(1/ K). 

We consider the properties of the second term of Equation ([6]) in the remainder of the 
paper, completing the proof of the theorem at the end of Section [31 □ 

Theorem 11.51 can be interpreted in the same way as Theorem 1 of [9], that 'each user 
is able to achieve the same rate that he would achieve if he had the channel to himself 
with no interferers, half the time'. 

The theorem is presented under the assumptions of uniform Pr and P^ with deterministic 
fading for the sake of simplicity of exposition - we believe that the main result will be 
robust to relaxation of these conditions. In Section 14.11 we consider whether the and 
Pt need necessarily be uniform. In Section 14.21 we introduce a variant of the model 
with a random fading term. Although the theorem is based on probabilistic arguments, 
in Section 14.31 we describe an associated algorithm which bounds the sum capacity of 
arbitrary Gaussian interference networks. 

1.6 Relation to previous work 

As reviewed in more detail by Jafar pj, recently progress has been made in several 
directions towards understanding the capacity of Gaussian interference networks. 

In problems concerning networks with a large number of nodes, work of Gupta and 
Kumar [6] uses techniques based on Voronoi tesselations to establish scaling laws (see 
also Xue and Kumar [17] for a review of the information theoretical techniques that can 
be applied to this problem). Under a similar model of dense random network placements, 
though using the same points as both transmitters and receivers, Ozgiir, Leveque and 
Tse [121 EES] use a hierarchical scheme, where nodes are successively assembled into groups 
of increasing size, each group collectively acting as a MIMO transmitter or receiver, and 
restricting to transmissions at a common rate. Theorems 3.1 and 3.2 of [13] show that 
for any e > there exists a constant c e and a fixed constant c\ such that 

c^ 1 " < Cx < Cl K\ogK. (8) 

These bounds are close to stating that Cs grows like K, but without the explicit constant 
that Jafar [9] and Theorem II .51 of this paper achieve. In this paper, we produce a version 
of the upper bound of Equation (jHJ) without the logarithmic factor and being explicit 
about the constant ci, although note that this result is proved under a model that differs 
from that of [T3] in the fact that we have a total of 2K nodes rather than K. 

An alternative approach to Gaussian interference networks is to consider the limit of the 
capacity as the SNR tends to infinity, with a fixed number of users. Cadambe and Jafar 
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[2] used interference alignment to deduce the limiting behaviour within o(log(SNR)). 
These techniques were extended by the same authors [3] to more general models in the 
presence of feedback and other effects. 

For small networks, the classical bounds due to Han and Kobayashi [7] for the two-user 
Gaussian interference network have recently been extended and refined. For example 
Etkin, Tse and Wang [3] have produced a characterization of capacity accurate to within 
one bit. These results were extended by Bresler, Parekh and Tse [Tf, using insights 
based on a deterministic channel which approximates the Gaussian channel with suffi- 
cient accuracy, to prove results for many-to-one and one-to-many Gaussian interference 
channels. 

We briefly mention an alternative proof of Theorem 5 of [9], which provides a faster 
rate of convergence than that achieved in Equation (CQ). We first review some facts 
from graph theory, concerning random bipartite graphs formed by two sets of vertices 
of size N, with edges present independently with probability 5. Erdos and Renyi [I] 
prove that the probability of a complete matching failing to exist tends to for any 
5 = S(K) = (log if + ck)/K, where ck — > oo. We recall the argument where S is fixed, 
so that we can be precise about the bounds, rather than just working asymptotically. 
As in, for example, Walkup [IB] , we say that a subset As C A of size k and a subset 
Bs C B of size N — k + 1 form a blocking pair of size k if no edge of the graph connects 
As to Bs- Equation (1) of [16] uses Konig's theorem to deduce that 

N 

P(no matching A to B) < ^ ^ P (.(As, Bs) blocking pair) (9) 

fe=l \A s \=k,\B s \=N-k+l 

By splitting the sum into terms where k < \/~N and k > \^N, so that (1 — p) k ( N ~ k+1 ) j s 
bounded by exp(— p(N + l)/2) and exp(— pN 3 ^ 2 /2) respectively, a bound of 

2v / iVW 2 ^exp(-p(iV + l)/2) + 2 2N exp(- P N 3 / 2 /2) 

can be obtained. We deduce that the probability of a complete matching failing to exist 
decays at an exponential rate in N. 

To prove Theorem 5 of [9], we divide the receiver-transmitter links into two groups A 
and B of size N = \_K/2\, and consider complete matchings on the bipartite graph 
between them. Each edge is present in the bipartite graph if the corresponding INR 
lies in a particular range (see Lemma |2~T1 below for details), which occurs independently 
with probability 5. For each pair that is successfully matched up, the corresponding 
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two-user channel becomes an e-bottleneck, and contributes < log(l + 2SNR) + e to the 
sum capacity. Hence, the high probability of a complete matching implies exponential 
decay of F(\C S /K - \ log(l + 2SNR)| > e), improving the 0(K~ 2 ) rate in Equation (JH). 

2 Technical lemmas 

We continue to work towards our proof of Theorem 11.51 by establishing some technical 
results. First in Section |2~TI we identify a condition under which the sum capacity of the 
two user interference channel can be bounded. Next, in Section I2~2l we show that the 
spatial separation condition of Definition [L2] holds under a variety of conditions. Further 
we show that spatial separation and the decaying condition Definition !!. 31 together imply 
that no SNR can be 'too large'. 

2.1 Bounds on two user channel 

First, we identify a condition on the values of SNR and INR under which the capacity 
of the two user interference channel can be bounded. The proof of the following result is 
based on Lemma 1 of [9], which gave the key definition of a 'bottleneck state', deducing 
that in the case SNRj = SNR.,- = INRjj = SNR, the sum capacity equals log(l + 2SNR). 
We reproduce the argument used there, to deduce a stability result that allows us to 
deduce when an e-bottleneck state occurs. 

Lemma 2.1 For any i,j , consider the two user inferference channel defined by 

Yi = exp(i<j) U ) \/ SNRjXj + exp(i<^) ■ s /lKR~ i X j + Z i: 
Yj = exp (^)0NR~X,+ &xp(i(j)jj) SNRj-Xj + Zj, 

where Zi and Zj are IID standard complex Gaussians, and <pij are independent U[0, 2tt] 
random variables independent of all other terms. 

If INK ji > SNRj then any reliable transmission rates satisfy 

R\i) + R[j] < log(l + INRjj + SNR,). (10) 

Proof We adapt the argument of Lemma 1 of [9]. That is, consider reliable transmission 
rates R[i], R[j]- Since the transmissions are reliable, then receiver i can determine X, 
with an arbitrarily low probability of error. 

Again, reliable transmission rates would remain reliable if receiver j was presented with 
Xi by a genie. In that case, it is easier for receiver i to determine Xj than it is for 
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receiver j to do so (since the weighting INRjj is larger than SNRj). However, we know 
that receiver j can recover Xj, since the rate R[j] is reliable, so we deduce that receiver 
% must be able to do the same. 

Since receiver % can determine X^ and Xj reliably, then these messages must have been 
transmitted at a sum rate lower than the sum capacity of a two-user multiple access 
channel, see for example [151 Equation (6.6)], which is log(l + SNRj + INRjj). □ 



2.2 Decay of tails 



Recall that the node positioning model given in Definition 11.11 involves independent 
positions of nodes sampled from identical distributions Pt and P^. We give examples of 
conditions under which the spatial separation property of Definition 11.21 holds. 



Lemma 2.2 



(i) If either Ft orF R has a density with respect to Lebesgue measure which is bounded 
above on T> then Pt and F R are spatially separated. 

(ii) IfF T and F R are supported on sets T and TZ that are physically separated, in that 



d* = inf {d(t, r) : t G T, r G IZ} > 0, 



then F T and F R are spatially separated. 



Proof (i) If Ft has a density with respect to Lebesgue measure which is bounded above 
by C, then for any ball B s (y) of radius s centred on y, the probability Pt(-B s (?/)) < 
CVt>s d , where Vd is the volume of a Euclidean ball of unit radius in M D . Hence 



V(d(T,R) < s) = J F T (B s (y))dF R (y) < CV D s D j dF R (y) = CV D s D 

and the result follows, taking C sep = CVd and D sep = D. The corresponding result for 
F R follows on exchanging P^ and Pt in the displayed equation above. 



If s < d*, then F(d(T, R) < s) = 0. If s > 4, then we have F(d{T, R) < s) < 1 < 
s/df, so that we can take C sep = 1/df and D sep = D. □ 



Next we show that combining the spatial separation condition of Definition 11.21 and the 
decaying condition Definition 11.31 gives us good control of the maximum of Su. The 
argument is similar to that given in Theorem 3.1 of [13J. 
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Lemma 2.3 Consider an IID network with spatially separated P^ and P#, with signals 
decaying at rate a. The probability that the maximum of the K random variables Su is 
large tends to zero: 

lim P ( max S u > al ° gK \ = (11) 

K^oo \l<i<-K" D sep J 

where D sep is the separation exponent from Definition \l. 6 2[ 

Proof We combine Definitions 11.21 and 11.31 Since all SV, have the same marginal 
distribution, it is enough to deduce that, for any u > 1, 

P(5« > u) = P(l/2 log(l + 2SNRj) > u) 
= P(SNRi > (exp(2u) - l)/2) 

< P(SNRi > exp(2M)/3) 

= ¥(f(d(T u R i ))>exp(2u)/3) 

< P(C dcc /(rf(T,, > exp(2u)/3) 

= P(rf(T 4 , J R l ) < (C dcc /3) 1/a exp(-2u/a)) 

< C scp (C dcc /3) D ^ a exp(-2uD sep /a)), (12) 

where C sep is the separation constant from Equation (jHJ) and Cd ec is the decay constant 
from Equation (jlj). The result follows from Equation ( fT2|) using the union bound: 

P | max & > < AT f S „ > < ^(CW3) W " 



i<i<^ D sep y y D scp j k 2 

which tends to zero as K — > oo. □ 



3 Proof of Theorem 11.51 



In this section we complete the proof of the upper bound in Theorem 11.51 We calculate 
bounds on sum capacity using a strategy suggested by the work of Jafar [S] , who proves 
that the presence of a large number of disjoint two-user channels, each close to being a 
'bottleneck state', allows good control of the channel capacity. 

First in Section 13.11 we partition space into regions B uv . Lemma 13.21 gives an upper 
bound on the sum capacity of a two-user channel made up of points in neighbouring 
regions. Further, Lemma 13.31 tells us that we can control the number of links in each 
region. In Section l3~2l we complete the argument by matching elements of box B u v with 
elements of -B u _ e , v +e- This allows us to control the overall sum capacity of the K users, 
and to complete the proof of Theorem 11.51 
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3.1 Spatial partitioning model 



We consider each transmitter-receiver pair as a point in the joint domain VxD. Each 
transmitter-receiver pair (Tj, Ri) G T> x T> can be placed in well-defined disjoint regions 
B u v , allowing us to control the performance of the corresponding link. We write x^ for 
the Ith. coordinate of the vector x = (x^\ . . . , x^). 

Definition 3.1 Given M, we partition the space M. 2D and hence the joint transmitter- 
receiver domain DxD by a regular grid of spacings 1/M. For each u G IP and v G IP , 
we define boxes labelled by their 'bottom-left' corner as 

«~ " (M = ^^<^^<^ *r-l}.M) 

We wnte 5 = {(u, v) : B uv f)(V x D) ^ 3} /or £/ie se£ of possible labels (u, v), and 
split T> x T> into orthants, indexed by vectors E(u, v) G {—1,0, 1} D , with co-ordinate 

if v ii) _ w (0 > ; 

£«(u,v) = { ifv®-u® = 0, 

if v {i) _ u (i) < o. 

We introduce two subsets of S that we will not attempt to match, according to the rule 
that B u v is matched with B u - e v+e , where e = E(u,\). 

(i) First, the 'spine' S spinc = {(u, v) : E®(u,v) = for some I}. 

(ii) Second, the 'edge' iS e d ge = {( u > v ) : (5 U , V IJ -B u _ e v+e ) G\ V x V}. That is, the 
regions which overlap the boundary ofVxV, or which are matched with a region 
that overlaps the boundary. 

The control of position obtained by matching links in two neighbouring regions converts 
into control of the values of SNR and INR, allowing the sum capacity of the two pairs 
to be bounded using Lemma 12.11 Figure [2] gives a schematic diagram of a pair of 
boxes which we attempt to match using the construction in Lemma 13.21 plotting both 
transmitter and receiver positions singly on T> rather than jointly onVxT>. 




Lemma 3.2 Suppose the receiver-transmitter pair (i?,,Tj) appears in region B u v and 
the receiver-transmitter pair (Rj, Tj) appears in region B u _ e v+e , where (u, v) G iS\5 S p ine 
and e = E(u, v), then any reliable rates for those two links satisfy 

R[i\+R\j]<log(l + 2f(d(u,v))), 

where d(u, v) is the minimum transmitter-receiver distance between (Ri,Ti) G B u v . 
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Figure 2: Schematic plot of matched boxes, where boxes are labelled according to their 
bottom-left corner. We consider the case where D = [0, l] 2 , and plot the receiver and 
transmitter positions on the same square, with e = -E^u, v) = (—1, 1). The key property 
is that we can observe that d(Tj,Rj) > d(Tj,Ri) > d(Ti,Ri). We prove this rigorously 
in order to prove Lemma 13.21 
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Proof The rates are only improved by being presented with the messages of all the other 
users, reducing the situation to that of Lemma 12.11 For each I such that co-ordinate 
e« = 1, by construction (a) Tf - Bf > 0, (b) Rf < u« < Rf, (c) > u« > Tf. 
Hence, by (b), 

(Tf - Rf) = (Tf - Rf) + (Rf - Rf) < (Tf - Rf), 

and by (c), 

(Tf - Rf) = (Tf - Rf) + (Tf - T? y ) > (Tf - Rf). 
Overall then, in the case = 1, 

< (T^ l) - Rf) < (Tf - Rf) < (Tf - Rf). 

A similar argument applies for each I with = — 1, with the order of the signs reversed. 

Overall, we deduce that d(Tj, Rj) > d(Tj, Ri) > d(Ti, R4) > d(u, v), or that f(d(u, v)) > 
SNRj > INRjj > SNRj, so that Lemma [2.11 applies, allowing us to bound 

R[i] + R\j] < log(l + INRjj + SNRj) < log (1 + 2f(d(u, v))) , 

as required. □ 

Each of the K transmitter-receiver links are placed in the regions B uv where (u, v) e S. 
Our node positioning rule implies that each link is placed in region B uv with some 
probability p U]V independently of the others. We write iV UjV for the total number of links 
placed in region B u v , noting that the marginal distribution of each iV U)V is Bin(K,p uv ). 

Lemma 3.3 We can bound the probability that any of the regions contain a significantly 
different number of links to that expected at random: 

¥ ( max |JV UV - ^Puvl > kA < K l ~ 2r] . 
V(u,v)es' ' ' / 

Proof A standard argument using the union bound and Chebyshev gives 

F(max s \N U:V -K Pu , v \>KA = P I |J {\N U , V - K Pu , v \ > K*}\ 

V(U ' V)6,S ' \(u,v)SS / 

< Yl P (l^u,v-^u,v| >tC>) 

(u,v)es 



(u,v)e5 



Var (N^ 
K 2r > 
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since Var (iV U)V ) = Kp u>y (l -Pu,v) < %,v, so that J2( u ,v)es Var ( N u,v) <K. □ 



3.2 Matching links 

We now complete the proof of Theorem 11.51 - recall that we consider uniform node 
distributions P-r and P^ona bounded domain T> with smooth boundary. 

Proof of Theorem 11.51 The total sum capacity Ce < Im + Jm, where Im is the 
contribution from matched pairs of links and Jm is the contribution from unmatched 
links. We will consider M growing as a power of K, but for now, it is enough to regard 
M as fixed. 

We pair up the remaining edges in S \ (5 sp ine U <S c dge), working orthant by orthant. In 
particular, the matching between u, v and u — e, v + e is one-to-one, and each region is 
counted at most once. 

For each e £ { — 1, 1} D , we can define the function Il e by Il e (w) = w ■ e. The key 
observation is that for each (u, v) ^ 5 S pme, if e = -E^u, v) then the inner product 
< n e (v - u) < n e ((v + e) - (u - e))n e (u - v) + 2D, so that (u - e, v + e) ^ S spinc . 
We can sort the regions B uv by value of n e (u — v), at each stage adding some (u, v) 
with the lowest value of II e (u — v) that has not yet been matched to the set <Sb dy We 
depict this matching in Figure [3 
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Figure 3: Partition of regions into <S> spine , 5 e d ge and iSbody This illustrates the case where 
V = [0, 1] and M = 6, and we plot the regions as squares on V x V = [0, l] 2 . We label 
regions in iS sp i ne by S, regions in iS c d gc by E and regions in <Sbod y by B, with an arrow 
into the region they are matched with. 
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This means that we can select a set of regions <Sbod y sucri that for u, v e <Sbod y , the ^u,v 
and B u _ e v+e between them cover S \ (5 S p; nc U <S e dge); that is 



U (£ u ,v u B, 




u,ve5 body 



For V with volume V and each (u, v) G <5>bod y the probabilities p UjV = £> u -e,v+e = 
1/(V a2 M 2jD ), since the two boxes do not intersect the boundary of D x T>. By assumption 
there are at least Kj (V 2 M 2D ) — K v links in each of the regions -B u ,v and -B u _ e , v +e, since 
Lemma [3.31 shows that the probability that this does not occur is < K 1 " 271 . We choose 
K/(V 2 M 2D ) — K 71 links randomly from each pair of sets and match them, with Lemma 
13.21 implying that each pair of matched links contributes at most log(l + 2/(d(u, v))) 
to the sum capacity. Overall, the total contribution to the sum capacity from all the 
matched links satisfies 



By the definition of Riemann integration, by picking M sufficiently large, the term 
hi/K < ~Elog(l + 2SNR) + e/2, for e arbitrarily small. 

Next we control Jm, the contribution from the unmatched links. Specifically, without 
loss of generality, if D is bounded with volume V, we can assume T> x T> C [0, L] 2D for 
some L. 

(i) We do not attempt to match some regions because they belong in S e d ge or iS sp i ne . 

(a) We assume that the boundary of T> is sufficiently smooth that there exists a 
finite A ('surface area') such that |iS e dgc| — AM 2D ~ X for all M. (For example 
if V = [0, 1] D then |5 edge | < 4DM 2D ~ l , since there are 2D co-ordinates that 
can take values or M — 1, and then M 2D ~ l values for the remaining co- 
ordinates) . 

(b) The number of regions in |iS sp i ne | < D(LM) 2D ~ l , since there are D co- 
ordinates which can agree, and at most LM possible values the remaining 
co-ordinates can take. 

Overall, there are at most (A+DL 20 ' 1 )^! 20 " 1 regions we do not attempt to match. 
Each region we do not attempt to match contains at most Kj (V 2 M 2D ) + K n links. 

(ii) We attempt to perform matching between at most (LM) 2D regions, with at most 
2K n unmatched links remaining from each. 




(14) 



u,veS body 
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H J m/K > e/2) < P max S u > 



In total we deduce there are (A + DL 2D - l )(K) \V 2 M) + K^M 2D ~ l ) + 2(LM) 2D K r > un- 
matched links. Writing (3 = 1/(3 (2D + 1)), and choosing M = i^ 3 ^( 1_r ?) ; and for example 
taking i] = 2/3, there are (K 1-13 ) unmatched links. 

The single user capacity bound (see for example [H3, Equation (6.4)]) tells us that reliable 
rates satisfy: 

R\i) < log(l + SNRj) < log(l + 2SNRj) = 2S U < 2 max^, (15) 

i 

so that Jm = C\K l ~P max!<j<^ Su for some c\. Hence overall, the probability 

eg ' 

i<i<k '~ " ~ 2ci 

which tends to zero by Equation ( iTTl) . □ 

Note that the upper bound in Equation ([H]) obtained by Ozgiir, Leveque, and Tse [T3] 
essentially has the extra factor of log K since the bound is only made up on the term 
Jm- It is precisely the matching argument that gives rise to the Im term which has 
reduced the order of the bound, as we take advantage of the extra randomness provided 
by placing transmitter and receiver nodes separately. 

Note that Equations ( TTTT) and (|T5l) together give probabilistic bounds on maxi<j<^- R[i}. 
Specifically, since they prove that max!<j<x R[i] = Op(logi^), they control the extent 
to which a large sum capacity can be achieved by a small number of links that operate 
at particularly high capacity. This suggests that in this case the sum capacity is not too 
unfair a measure of network performance. 



4 Future work and extensions 

We briefly comment on some extensions of Theorem 11.51 to more general models. 



4.1 Non- uniform node distributions 

We would like to consider more general distributions P# and Ft, rather than simply 
assuming that these distributions are uniform. Indeed, we might wish to extend to a 
situation where P^t have a joint distribution from which we sample independently to 
find receiver and transmitter distributions. 

The main issue that arises is whether we can quantize the joint distribution into regions 
_B U)V which can be uniquely paired off with other regions B u iy, with the paired regions 
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having equal probability, and with bounds on the sum of reliable rates being possible 
according to results such as Lemma 13.21 

One case that certainly works is that where the co-ordinates of P# and Ft are indepen- 
dent. In this case we can extend Equation (fl3l) to obtain 



where qn,i(%) is the xth quantile of the distribution of the /th component of the re- 
ceivers, and qr,i(x) is the xth quantile of the distribution of the /th component of the 
transmitters. 

The proof of Theorem 11.51 follows exactly as before in this case. 
4.2 Random fading model 

We briefly remark on an adaption of Definition 11.41 that would have the same indepen- 
dence structure, while introducing random fading into the model. That is, we could 

set 



where the My are IID random variables with a density. Under our node placement 
model, the SNRj will again be IID, making this model tractable in much the same way. 
In particular we can extend the tail behaviour bounds given in Lemma [2.31 to this case, 
adjusting the constant to take account of the random fading term. 

Lemma 4.1 Consider an IID network with spatially separated Ft and Fr, with sig- 
nals decaying at rate a. If the fading random variables M y - have finite mean then the 
probability that the maximum of the K random variables Sa is large tends to zero: 



B 



u,v 



{(x,y) : q n ,i(u/M) < x® < qiZ)l ((u + 1)/M), 

<tr,i(y/M) < y w < qr tl ((v + 1)/M), for all /} , 



INRjj = Mijf(d(ti,rj)), 




Proof An equivalent of Equation ( fl2l) holds, since as before, for any u > 1, 



F(S U > u) < P(M ii /(d(T ijJ R i )) > exp(2«)/3) 




< P(M« > exp(w)/V3) + F{C dec /d{T h Ri) a > exp(u)/V3) 
= VSEM U exp(-u) + F(d(T h R t ) < (C dec V3) 1/a exp(-u/a)) 

< y/SKMa exp(- M ) + C scp (C dcc V3) D - /a exp(-uD scp /a)). 



18 



so the result follows exactly as before, using the union bound. 



□ 



The key to proving convergence in probability of C-^/K is to show that matching is 
possible between elements of B u v and -B u -e,v+ e) as before. In the case of deterministic 
fading, any links (i?j,T,) G B uv and (Rj,Tj) G B u _ ev+e could be matched, since the 
proof of Lemma 13.21 showed that in this case 

d{Tj, R 3 ) > d(Tj, Ri) > d(T h Ri). (16) 

In the case of random fading, this is not enough to control the relevant values of INR. 
However, Lemma 14.21 shows that we can match a high proportion of links, by looking 
for (i, j) such that 

Mjj < Mji < M u , (17) 

which can be combined with Equation ( TTBl to deduce that SNRj < INR^j < SNRj, so 
that again the sum capacity of the relevant two user channel < log(l + 2SNRj). Note 
that it is enough for our purposes to consider the case of uniform My with densities, 
since only the ordering between random variables matters in Equation f|T7|) . We give a 
technical lemma that will imply the control that we require. 



Lemma 4.2 Consider a bipartite graph with n vertices in each part which we refer to 
as A= (Ai, . . . , A n ) and B = (Bi, . . . , B n ) respectively. All the vertices are labelled with 
independent U[0, 1] random variables, with Ai labelled by Ui and Bj labelled by Vj. The 
bipartite graph has an edge from Ai to Bj iff Ui < Wij < Vj, where Wij are U[0,1], 
independent of (U, V) and each other. 

For any 7 > 2/3, there exists a matching of all but 0(n 7 ) vertices, with probability 
> 1 - 5exp(-2n 2 ^ 1 ). 



Proof We throw away the 3n 7 vertices with the biggest values of Ui and the 3n 7 
vertices with the lowest values of Vj. This leaves new sets A and B each of size N = 
N(n) = n — 3n J . We will show that there exists a matching between A and B with high 
probability, again using Equation @ and controlling the probability of blocking pairs. 
We condition on the event 



sup 

t 



#0 : Ui < t} 



t 



11 



< n 1 



>}U{sup 



#0 : Vi < t} 



t 



n 



< n 



7-1 



since Massart's form of the Dvoretzky-Kiefer-Wolfowitz theorem [TU] tells us that this 
does not take place with probability < 4 exp(— 2n 27_1 ). 



Conditional on the event ( jl8l) . for any k we know there are more than n — k — 2n 7 values 
of Ui which are less than 1 — k/n — n 1 ^ 1 . Equivalently, there are fewer than k + 2n 7 
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values of Ui larger than 1 — k/n — n 1 ^ 1 . Since we throw away the largest 3n 7 values 
of A, any subset of As C A of size k has at least n 7 vertices with Z7 f values less than 
1 — k/n — n 1 ^ 1 ('small vertices'). 

By a similar argument, any subset Bs C 23 of size iV — + 1 has at least n 1 vertices with 
Vj values greater than 1 — k/n + n 1+1 ('large vertices'). See Figure 0] for a depiction of 
these events. 



U- 



xxxxxxxxxxxxxxxxxxxxxxxiaaa )qq< xx 



x^xx^xxxxxxxxxxxxxxxxxxxxxxx 



1 — k/n 



n 



7-1 



1 — k/n + n 



7-1 



1 — k/n 



Figure 4: Position of vertices in subsets As and Bs forced by throwing away 3n 7 largest 
values of U and 3n 7 smallest values of V . 



There is an edge between each of these small vertices in As and large vertices in Bs 
independently with probability at least 2n J ~ 1 . Hence, the probability that a particular 
As and Bs form a blocking pair is less than (l-2n<~ 1 ) n 7 < exp(— 2n 37 1 ). Substituting 
in Equation (Q the probability of no matching is 

^ E (I) ( N f k + x) ex P (-2n 37 -) = exp(-2^-i) < 2 « ex p(-2n^ 1 ), 

and the result follows since 37 — 1 > 1 , combining with the probability of (1181) failing to 
occur. □ 

The remainder of the proof of Theorem 11.51 carries over as before. We need only alter 
Section I3~2l and this can be done the increase in numbers of unmatched vertices remains 
sublinear in K. 
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4.3 Constructive algorithm 



Although Theorem 11.51 only gives a result concerning average performance of large net- 
works, it does suggest some techniques that can be used to approximate the sum capacity 
of any particular Gaussian interference network. If the network is created via a spatial 
model, then we can attempt to match cross-links into e-bottleneck channels using the 
constraints on spatial position described in the proof of Theorem 11.51 deducing bounds 
clS cL result. 

However, even if we are only presented with the values of SNRj and INRjj, it may be 
possible to find bounds on sum capacity using the insights given by Lemma 12.11 Using 
the interference alignment scheme of [TT] , we know that a lower bound on Cs is given 
by Ei>g(l + 2SNR,). 

One possible algorithm to find an upper bound works as follows: 

(i) Sort the indices by value of SNRj, and for some M, partition the transmitter- 
receiver links into 2M categories B r of approximately equal size K/(2M). 

(ii) For each M, we attempt to match links between B 2 m-i and B 2 m, considering the 
bipartite graph between them. 

(a) We add an edge to the bipartite graph between j e B 2 m-i and i G B 2M if 
SNRy < INRjj < SNRj. 

(b) We look for a maximal matching on the bipartite graph, using (for example) 
the Hopcroft-Karp algorithm [8j, which has complexity y/VE, where V is the 
number of vertices and E the number of edges. 

(iii) By Lemma 12.11 each edge (j, i) in each maximal matching contributes log(l + 
INRjj + SNRj) as an upper bound on the sum capacity, and each unmatched 
vertex i simply contributes the single user upper bound of log(l + SNRj) 

By varying the size of M, this algorithm will find a range of upper bounds, of which 
we can choose the tightest. We want M large enough that categories B r each contain 
a narrow range of SNR values, but M small enough that there are plenty of points in 
each range B r to ensure a large maximal matching. 

5 Conclusions 

In this paper, we have deduced sharp bounds for the sum capacity of a Gaussian inter- 
ference network. Our main contribution comes through the upper bound, which uses 
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arguments based on controlling the position of pairs of vertices to match them into bot- 
tleneck states. Although our main result is proved under the assumption of deterministic 
fading, with signal strength decaying as a function of distance, in Section fl~2l we describe 
an extension to a model with random fading. 
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